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c : 

Abstract 

o : 

£f~) , Recent developments in the field of automatic computation of Feynman diagrams using 

asymptotic expansions are reviewed. The hadronic decay rate of the Z-boson is taken as an 
' ' , example for their physical application. 

^ ■ 1 Introduction 

: 

The final round of the analysis of data taken at LEP in the runs between 1990 and 1995 at 
energies around the Z-peak is going to be completed 0. About 16 millions of Z-bosons have 
been produced and the resulting accuracy of, e.g., the Z-boson mass and its decay width is 
impressive (for the latest experimental values of these quantities see Q). In view of this precision 
there have been huge efforts in calculating higher order corrections to these observables. This 
talk is supposed to be not so much a review of theoretical and experimental results but is rather 
i-G ■ concerned with some of the available tools to perform such calculations. 

X ; 2 QCD corrections and asymptotic expansions 



The decay rate of the Z-boson into quarks constitutes an illuminating example of how to use 
asymptotic expansions of Feynman diagrams to simplify calculations. In what follows we will 
only be concerned with the vertex corrections and it will always be assumed that they are com- 
puted via the optical theorem by calculating the Z-boson self energy and taking the imaginary 
part. If one considers, e.g., only QCD-corrections, it certainly is a reasonable lowest order ap- 
proximation to neglect all quark masses. Only one dimensional quantity is left, the Z-boson 
mass, and one arrives at massless propagator diagrams for which the so-called integration-by- 
parts algorithm is available and has been explicitely worked out |2| and implemented in a FORM 
H package called MINCER Q up to three loops. This immediately gives the answer for the 
Z-boson decay rate into massless quarks up to 0(a 2 s ) ||. Making additionally use of infrared 
re-arrangement, it is even possible to extend the result to 0(a^) || 0- 

As a second step one may want to take effects induced by the quark masses into acccount. 
Then, however, one is faced with two-scale Feynman integrals. Although the full 0(a s )-result is 
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Figure 1: Subdiagrams contributing to the Large Momentum Procedure of a particular diagram. 



known j|, [| 10] — the vector part even for a very long time — , at 0(a 2 ) only certain subclasses 



of diagrams have been computed analytically (e.g. IT] 

For the remaining part one is forced to consider some approximation procedure to obtain, 
e.g., an expansion in m 2 /M z , with m q the mass of the quark in the final state. The recipe 
for the efficient computation of such asymptotic expansions has been worked out in a series of 
publications and shall not be repeated here (for a review see [[12;1). The essence is to expand 
the integrands of certain subgraphs of the initial diagrams, leading to a complete factorization 
of "small" and "large" quantities. Two special cases may be distinguished: When only masses 
appear as large quantities, the technique is called the Hard Mass Procedure. In contrast, the 
Large Momentum Procedure deals with the case of only large momenta. 

Concerning the quark mass effects in the QCD-corrections to the Z-boson decay rate, it is 
clear that here the appropriate procedure is the second one. There is only one large (q 2 = Mf ) 
and one small (m 2 ) mass scale in this problem, so the factorization mentioned above means that 
products of only single scale diagrams are produced: massless propagator diagrams (m = and 
g / 0) and massive tadpoles (to ^ and q = 0). We already mentioned the integration- by- 
parts algorithm to compute the former ones. The underlying principle may also be applied to 
the latter ones, as has been done in [jjUfl , again up to three loops. The implementation of this 
procedure has been performed, for example, in a FORM package named MATAD [14]. 



There are, of course, certain limitations of this approach both from the technical and the 
analytical point of view. The former one is connected with the realization of the prescriptions 
provided by the asymptotic expansions. To three-loop order it becomes a non-trivial task to find 
and properly expand the contributing subdiagrams. For one particular diagram, all subgraphs 
contributing to the Large Momentum Procedure are shown in Fig. [|. The solution of this 
problem is to use the algorithmic nature of the prescriptions and to pass their evaluation to a 
computer. As far as the Large Momentum Procedure for two-point functions is concerned this 



has been performed in a PERL program called LMP [15]. Another program, doing also the Hard 
Mass Procedure, will be mentioned in Section ||. Therefore, in this sense this technical limitation 
no longer exists. 

The second limitation of asymptotic expansions is more severe. It is clear that a cut series in 
general contains less information than the full result. In our case of the Z-boson decay, let us, 
for example, consider the hypothetical case 2m q < Mz < ^m q . The production of four fermions 
is then kinematically forbidden. A small m^-expansion, however, is uncapable of discriminating 
among any of the cases Mz > m q , Mz > 2m q or Mz > 4:m q , which is why at first sight 
one cannot expect to obtain reasonable results below the four-quark threshold. However, the 
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Figure 2: Non-abelian contribution to i?had as an example for the application of the Large 
Momentum Procedure |jq| , plotted over x = 2m/y / i. Narrow dots: semi-analytical result from 
17], other lines: power corrections up to (m?/s) n , n = 0, . . . ,6. 



situation is not so bad as it seems. The reason is that in general the four-particle channel has a 
very smooth threshold behaviour due to phase space suppression. In a small-m g expansion this 
smoothness is carried over to energies below the four-particle threshold preserving its validity 
in this region [|l6| . 

As an example the contribution to the hadronic .R-ratio proportional to the color factor 
CaCf as a function of x = 2m (J /y / s, where s is the cms-energy, is shown in Fig. ||[ The first 
threshold is at yfs = 2m q {x = 1), the second one at y's = Am g {x = 1/2). Convergence, 
however, seems to be warranted at least up to x = 0.8. 

Anyway, for the Z-boson decay one is above any possible four-quark threshold such that 
convergence of the expansion to the correct answer is guaranteed. Not only this: The quadratic 
and quartic mass corrections, which already provide a very good approximation in this case, 



can also be obtained without really using asymptotic expansions, even up to O(a^) |L8], 15]. As 
far as the Large Momentum Procedure is concerned, there certainly are more important fields 
of application for it (see, e.g., [JH], |(J). On the other hand, for the Hard Mass Procedure, 
important contributions to Z-decay have recently been computed where it really proves to be a 
very useful tool, as we will see in Section ||. 

3 Electroweak and mixed QED/QCD-corrections 

Let us turn to the electroweak radiative corrections to the Z-boson decay. As far as one only 
considers QED, the similarity to QCD allows not only to evaluate all pure QED-corrections, 



but also QED/QCD-corrections of order aa s , a 2 a s and ac? s [21] from the knowledge of the 
diagram-wise results to 0{a 2 s ) resp. 0(a^) by altering the color factors. 

The situation is different when allowing also for Z- and VK-boson exchange between the 
produced quarks, because of their non-negligible masses. As far as the W-boson is concerned, 
another phenomenon occurs, namely the appearence of the isospin partners of the produced 
quarks as virtual particles in the loops. In contrast to the decay into u-, d-, s- or c-quarks, 
where this does not really produce a difference in comparison to virtual Z-exchange because all 
of them may be considered as massless, the decay into 6-quarks is somewhat exceptional because 
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the 6-quark is the isospin partner of the i-quark. Neglecting the i-quark mass certainly is not 
a good approximation, nor can it be set infinite as one knows, for example, from the radiative 
corrections to the p-parameter that appear to be proportional to m 2 |2^ ]. Nevertheless, a full 
result for T(Z — * qq) to 0(a) for the decay into q = u,d,s,c |23|, |||] as well as for the one into b 
p5| , p4| is available. The leading m^-behavior for the latter is quadratic like for the p- parameter. 
To 0(a 2 ) the i-quark enters also the calculation for the decay into u,d,s,c. The leading mf— 
and the subleading m^-terms to this order are known [26, 27]. For the decay into b- quarks only 
the leading m^-terms are available p3|. 



4 Mixed electroweak/QCD-corrections 



A part of the mixed aa s -corrections, namely those induced by virtual gluon and photon ex- 
change, has already been mentioned in Section [|. The present section will be concerned with 
the case of W- or Z-, accompanied by an additional gluon-exchange, i.e., in a sense, QCD- 
corrections to the 0(a)-results described in Section ||. Again it is natural to distinguish between 
the decay into u,d,s,c and into b. The former case was evaluated in [28], and it is instructive 
to dwell a bit on the technique which was used in this work. 

Again the rate was determined by computing the Z-boson self energy up to the order con- 
sidered and taking the imaginary part of the result. In the case of virtual Z-boson exchange, 
one arrives at three-loop on-shell integrals, for Vl^-exchange one gets two-scale diagrams. The 
idea in |2g] was to use the Hard Mass Procedure described in Section |2| to expand the diagrams 
in M|/M 2 , where M is the mass of the virtual gauge boson, and take the limit M — > Mz resp. 
M — ► My/ in the final result. Since convergence at these points turned out to be quite slow, 
a part of the diagrams was also expanded in M 2 /M§. In this way it was possible to obtain a 
reasonable approximation to the full result. 

In the case of Z — > bb one faces the problem of an additional mass scale, the f-quark mass. 
Using the Hard Mass Procedure for m 2 M Z ,M^, one may factor out the m^-dependence. 
However, for a part of the diagrams one still is left with two-scale and even three-scale integrals 
involving M§ and and £iyM^, where £w is the electroweak gauge parameter which we want 
to keep. Although they appear to be only one-loop integrals, their exact evaluation up to 0(e) 



produces inconvenient results. Instead, the results of 2£] were obtained by applying the Hard 
Mass Procedure to these kinds of diagrams once more, this time using £iyM^, 3> M§. This 
seemingly unrealistic choice of scales becomes justified by recalling the discussion of Section ||: 
It is not possible for an expansion to distinguish the inequality 3> M§ from 4M^ S> M§ or 
{m% + M\y) 2 M§, the latter ones being perfectly alright. The only matter is to perform the 
expansion on the appropriate side of all thresholds, and here one is concerned with thresholds 
at 2Mw and at mt + My/. Therefore, the choice 3> M§ is to be understood purely in this 
technical sense. Graphically this continued expansion looks as follows: 




+ 



where only those terms are displayed which are relevant in the discussion above and all others 
contributing to the Hard Mass Procedure are merged into the ellipse. The thick plain line is the 
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top quark, the thick wavy one a Goldstone boson with mass squared £ypMjjp, for example. The 
thin plain lines are 6-quarks, the inner thin wavy lines are W-bosons, the outer ones Z-bosons. 
The spring-line is a gluon. The mass hierarchy is assumed to be m 2 S> £yyMyy S> 3> Mf . 
The freedom in choosing the magnitude of provides a welcome check of the routines and the 
results. 

The outcome of this procedure is a nested series: The coefficients of the M\y/mt -expansion 
are in turn series in Mz/M\y. Note that in contrast to the decay into u,d,s,c there is no 
threshold at M\y which makes an additional expansion in M\y/Mz unnecessary. 

In view of this calculation the procedure of successive application of the Hard Mass Procedure 
resp. the Large Momentum Procedure has been implemented in a Fortran 90 program named 



EXP [30]. Therefore, the computation of a three-loop two-point function can now be done fully 
automatically given some arbitrary hierarchy of mass scales. Even more, the link to the Feynman 
diagram generator QGRAF in a common environment called GEFICOM j3^] allows to obtain 
the result of a whole physical process without any human interference except for specification 
of the process and final renormalization. 

Finally, let us present the result for the VK-induced corrections to the Z-decay rate 5T W {Z — * 
bb) in the form of the renormalization scheme independent difference to the decay rate into dd. 
Inserting the on-shell top mass nit = 175 GeV, the Z-mass Mz = 91.91 GeV and sin 2 6y/ = 0.223 
gives 

8T W (Z bb) - 5T W {Z -► dd) = T ^ — -( - 0.50 

sin z Ow TT I 

+ (0.71 - 0.48) + (0.08 - 0.29) + (-0.01 - 0.07) + (-0.007 - 0.006) 
+ — 1.16 + (1.21 - 0.49) + (0.30 - 0.65) + (0.02 - 0.21 + 0.01) 

TT 



+ (-0.01 -0.04 + 0.004) 



^(- 0.50 - 0.07+^ 

Sin 0y/ TT I TT 



1.16 + 0.13 



(1) 



where the factor T° a /(tt sin 2 9w) with T° = Mz a/ (4 sin 2 0yy cos 2 6w) has been pulled out for 
convenience. The numbers after the first equality sign correspond to successively increasing 
orders in 1/m 2 , where the brackets collect the corresponding constant, logmj and, if present, 
log 2 mj-terms. The numbers after the second equality sign represent the leading m 2 -term and 
the sum of the subleading ones. The 0(a) and 0(aa s )-results are displayed separately. Com- 
parison of this expansion of the one-loop terms to the exact result of |24j] shows agreement up 
to 0.01% which gives quite some confidence in the aa s -contribution. One can see that although 
the m 2 -, m®- and m® log mt-terms are of the same order of magnitude, the final result is sur- 
prisingly well represented by the m 2 -term, since the subleading terms largely cancel among each 
other. 
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